TTPOO-10 - 



hep-ph/000523^ 



1 



Using differential equations to compute two-loop box integrals* 
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The calculation of exclusive observables beyond the one-loop level requires elaborate techniques for the compu- 
tation of multi-leg two-loop integrals. We discuss how the large number of different integrals appearing in actual 
two-loop calculations can be reduced to a small number of master integrals. An efficient method to compute these 
master integrals is to derive and solve differential equations in the external invariants for them. As an application 
of the differential equation method, we compute the C(e)-term of a particular combination of on-shell massless 
planar double box integrals, which appears in the tensor reduction of 2 — » 2 scattering amplitudes at two loops. 



1. Introduction 

Perturbative corrections to many inclusive 
quantities have been computed to the two- and 
three-loop level in past years. From the tech- 
nical point of view, these inclusive calculations 
correspond to the computation of multi-loop two- 
point functions, for which many elaborate calcu- 
lational tools have been developed. In contrast, 
corrections to exclusive observables, such as jet 
production rates, could up to now only be com- 
puted at the one-loop level. These calculations 
demand the computation of multi-leg amplitudes 
to the required number of loops, which beyond 
the one-loop level turn out to be a calculational 
challenge obstructing further progress. Despite 
considerable progress made in recent times, many 
of the two-loop integrals relevant for the calcu- 
lation of jet observables beyond next-to-leading 
order are still unknown. One particular class of 
yet unknown integrals appearing in the two-loop 
corrections to three jet production in electron- 
positron collisions, to two-plus-one jet production 
in electron-proton collisions and to vector boson 
plus jet production in proton-proton collisions are 
two-loop four-point functions with massless inter- 
nal propagators and one external leg off-shell. 

We elaborate on in this talk several tech- 
niques to compute multi-leg amplitudes beyond 



one loop. We demonstrate how integration-by- 
parts identities (already known to be a very valu- 
able tool in inclusive calculations) and identities 
following from Lorentz-invariance (which are non- 
trivial only for integrals depending on at least two 
independent external momenta) can be used to 
reduce the large number of different integrals ap- 
pearing in an actual calculation to a small number 
of master integrals. This reduction can be car- 
ried out mechanically (by means of a small chain 
of computer programs), without explicit reference 
to the actual structure of the integrals under con- 
sideration and can also be used for the reduction 
of tensor integrals beyond one loop. 

The master integrals themselves, however, can 
not be computed from these identities. We derive 
differential equations in the external momenta for 
them. Solving these differential equations, it is 
possible to compute the master integrals without 
explicitly carrying out any loop integration, so 
that this technique appears to be a valuable al- 
ternative to conventional approaches for the com- 
putation of multi-loop integrals. 

We demonstrate the application of these tools 
on several examples. 

2. Reduction to master integrals 

Any scalar massless two-loop integral can be 
brought into the form 



*Presented at Loops and Legs in Quantum Field Theory, 
April 2000, Bastei, Germany 



/(pi, . . . ,p n ) = 



2 



d d k d d l 



{2n) d (2ir) d D" 11 . . . D\ 



~°1 



■ °q 



(1) 



where the Di are massless scalar propagators, 
depending on k, I and the external momenta 
pi , . . . , p n while Si are scalar products of a loop 
momentum with an external momentum or of the 
two loop momenta. The topology (interconnec- 
tion of propagators and external momenta) of the 
integral is uniquely determined by specifying the 
set (Di, . . . , D t ) of t different propagators in the 
graph. The integral itself is then specified by the 
powers TOj of all propagators and by the selec- 
tion (Si, . . . , S q ) of scalar products and their pow- 
ers (m, . . . , n q ) (all the rrii are positive integers 
greater or equal to 1, while the rii are greater or 
equal to 0). Integrals of the same topology with 
the same dimension r = Y] ■ rrii of the denomina- 
tor and same total number s = J^. Ui of scalar 
products are denoted as a class of integrals 7 t .,. jS . 
The integration measure and scalar products ap- 
pearing the above expression are in Minkowskian 
space, with the usual causal prescription for all 
propagators. The loop integrations are carried 
out for arbitrary space-time dimension d, which 
acts as a regulator for divergencies appearing due 
to the ultraviolet or infrared behaviour of the in- 
tegrand (dimensional regularisation, 

The number N{It,r,s] of the integrals grows 
quickly as r, s increase, but the integrals are 
related among themselves by various identities. 
One class of identities follows from the fact that 
the integral over the total derivative with respect 
to any loop momentum vanishes in dimensional 
regularisation 
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where J is any combination of propagators, scalar 
products and loop momentum vectors. J can be a 
vector or tensor of any rank. The resulting iden- 
tities HH are called integration- by-parts (IBP) 
identities. 

In addition to the IBP identities, one can also 
exploit the fact that all integrals under consider- 
ation are Lorentz scalars (or, perhaps more pre- 
cisely, "d-rotational" scalars) , which are invari- 
ant under a Lorentz (or <i-rotational) transforma- 
tion of the external momenta 101. These Lorentz 



invariance (LI) identities are obtained from: 
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In the case of two-loop four-point functions, 
one has a total of 13 equations (10 IBP + 3 LI) 
for each integrand corresponding to an integral 
of class It, r ,s, relating integrals of the same topol- 
ogy with up to s + 1 scalar products and r + 1 
denominators, plus integrals of simpler topologies 
(i.e. with a smaller number of different denomina- 
tors). The 13 identities obtained starting from an 
integral It,r,s do contain integrals of the following 
types: 

• It,r,s'- the integral itself. 

• It—i,r,s- simpler topology. 

• 7t,r+i, S ) ^t,r+i,a+i : same topology, more 
complicated than h,r,s- 

• 7t,r-i,s)it,r-i,«-i: same topology, simpler 
than It, r .s- 

Quite in general, single identities of the above 
kind can be used to obtain the reduction of 
h,r+i,s+i or 7t )r +i, a integrals in terms of It,r,s and 
simpler integrals - rather than to get information 
on the h,r,s themselves. 

If one considers the set of all the identities 
obtained starting from the integrand of all the 
N(It, r ,s) integrals of class It, r ,s, one obtains 
(^Vibp + N-L,i)N(It >r , s ) identities which contain 
N(I t .r+i,s+i) + N(I t ,r+i, s ) integrals of more com- 
plicated structure. It was first noticed by S. La- 
porta §] that with increasing r and s the number 
of identities grows faster than the number of new 
unknown integrals. As a consequence, if for a 
given i-topology one considers the set of all the 
possible equations obtained by considering all the 
integrands up to certain values r* , s* of r, s, for 
large enough r*,s* the resulting system of equa- 
tions, apparently overconstrained, can be used for 
expressing the more complicated integrals, with 
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greater values of r, s in terms of simpler ones, with 
smaller values of r, s. An automatic procedure to 
preform this reduction by means of computer al- 
gebra is discussed in more detail in || . 

For any given four-point two-loop topology, 
this procedure can result either in a reduction 
towards a small number (typically one or two) 
of integrals of the topology under consideration 
and integrals of simpler topology (less different 
denominators), or even in a complete reduction of 
all integrals of the topology under consideration 
towards integrals with simpler topology. Left- 
over integrals of the topology under considera- 
tion are called irreducible master integrals or just 
master integrals. 

3. Differential equations for master inte- 
grals 

The IBP and LI identities allow to express in- 
tegrals of the form (|l|) as a linear combination 
of a few master integrals, i.e. integrals which are 
not further reducible, but have to be computed 
by some different method. 

At present, the complete set of master inte- 
grals for massless on-shell two-loop box integrals 
is known analytically Jf| — |lQ[] up to finite terms in 
e = (4— d)/2. For massless two-loop box integrals 
with one off-shell leg, several topologies are yet 
to be computed analytically. A purely numerical 
approach for computing these integrals order by 
order in e has recently been proposed by Binoth 
and Heinrich in 

A method for the analytic computation of mas- 
ter integrals avoiding the explicit integration over 
the loop momenta is to derive differential equa- 
tions in internal propagator masses or in ex- 
ternal momenta for the master integral, and to 
solve these with appropriate boundary condi- 
tions. This method has first been suggested by 
Kotikov to relate loop integrals with inter- 
nal masses to massless loop integrals. It has been 
elaborated in detail and generalised to differen- 
tial equations in external momenta in jlj|; first 
applications were presented in [Q. In the case 
of four-point functions with one external off-shell 
leg and no internal masses, one has three inde- 
pendent invariants, resulting in three differential 



equations. 

The derivatives in the invariants Sj j — (pi+pj) 2 
can be expressed by derivatives in the external 
momenta: 
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It is evident that acting with the right hand 
sides of (|J) on a master integral It,t,o will, after 
interchange of derivative and integration, yield a 
combination of integrals of the same type as ap- 
pearing in the IBP and LI identities for It,t,o, 
including integrals of type It,t+i,i an d h,t+i,o- 
Consequently, the scalar derivatives (on left hand 
side of (f|)) of I t .t,o can be expressed by a lin- 
ear combination of integrals up to It,t+i,i and 
It,t+i,o- These can all be reduced (for topologies 
containing only one master integral) to It,t,o and 
to integrals of simpler topology by applying the 
IBP and LI identities. As a result, we obtain for 
the master integral It,t,o an inhomogeneous linear 
first order differential equation in each invariant. 
For topologies with more than one master inte- 
gral, one finds a coupled system of first order dif- 
ferential equations. The inhomogeneous term in 
these differential equations contains only topolo- 
gies simpler than It,t,Q, which are considered to 
be known if working in a bottom-up approach. 

The master integral It,t,o is obtained by match- 
ing the general solution of its differential equation 
to an appropriate boundary condition. Quite in 
general, finding a boundary condition is a simpler 
problem than evaluating the whole integral, since 
it depends on a smaller number of kinematical 
variables. In some cases, the boundary condi- 
tion can even be determined from the differential 
equation itself. 

In writing down the differential equations for 
all master integrals appearing in the reduction of 
two-loop four-point functions with up to one off- 
shell leg, one finds that all boundary conditions 
can be related to the following four one-scale two- 
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loop integrals: 



following combination of master integrals arises: 
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These integrals fulfill homogeneous differential 
equations in s; their normalisation can there- 
fore not be calculated in the differential equation 
method. They can however be evaluated explic- 
itly using Feynman parameters p"5[ . 

4. Applications 

Using the differential equation method, we 
computed in Q all master integrals with up to 
t = 5 different denominators that can appear 
in the reduction of two-loop four-point functions 
with one off-shell leg. 

Differential equations were also used in the 
computation of two-loop on-shell master integrals 
with planar Q and crossed |10| topologies. These 
topologies contain two master integrals each. The 
differential equations were applied only to obtain 
a relation among both master integrals, while one 
of the master integrals was calculated by Smirnov 
for the planar case and Tausk for the non- 
planar case || using a different method. 

At this workshop, Nigel Glover reported first 
progress towards the computation of massless 
on-shell 2^2 scattering amplitudes at two 
loops |f6| . In the reduction of planar amplitudes, 
it turned out that it is not sufficient to know 
the planar double box integrals up to their finite 
terms in e. In the tensor reduction procedure, the 
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where x = t/s. The functions K\(x, e) and 
K 2 {x,e) were computed only up to finite terms 
in e by Smirnov and Veretin in |9) . The structure 
of the up to now unknown C(e)-term was sub- 
ject of debate at this workshop. We did therefore 
decide after the workshop to attempt its compu- 
tation using the differential equation method. 

Starting from the differential equations for the 
massless double box integrals with one off-shell 
leg, which were obtained using the algorithms de- 
scribed in Section ||, we can derive the differen- 
tial equations for the above on-shell integrals by 
constraining s 12 + s 13 + s 23 — (see also [p~p|| ) . 
Using s and x as independent variables, we ob- 
tain two homogeneous differential equations in s 
(corresponding to a rescaling relation, Q) and 
two coupled inhomogeneous equations in x, which 
can be employed to compute the two master in- 
tegrals. Using the fact that the master integrals 
and their derivatives are regular in t at t — — s, 
the boundary condition at x = —I is inferred from 
the differential equations. This condition is how- 
ever not yet sufficient to determine the bound- 
ary conditions for both master integrals, since it 
only fixes one combination of them. It does how- 
ever turn out that it is sufficient to match the 
non-logarithmic terms in Ki(x, e), obtained by 
Smirnov in Q|, up to finite order in e to find the 
boundary condition for K\(x, e) + K 2 (x, e) up to 
0(e). The differential equations are then solved 
order by order in e by expressing the unknown 
as sum of harmonic polylogarithms |l7j . Requir- 
ing the differential equation to be fulfilled and 
the boundary conditions to be matched, the co- 
efficients of the harmonic polylogarithms in the 
ansatz can all be determined. 

Using the same normalisation factor as ||] , the 
coefficient of the C(e)-term of the above combi- 
nation of double box integrals reads: 
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5. Outlook 



We have demonstrated how techniques devel- 
oped for multi-loop calculation of two-point func- 
tions can be extended towards integrals with a 
larger number of external legs. In particular, we 
have shown that the use of differential equations 
in external invariants can be used as a powerful 
method to compute master integrals without car- 
rying out explicit loop integrations. As a first 
example of the application of these tools in prac- 



tice, we computed some up to now unknown two- 
loop four-point functions, relevant for jet calcu- 
lus beyond the next-to-leading order. The most 
important potential application of these tools is 
the yet outstanding derivation of two-loop virtual 
corrections to exclusive quantities, such as jet ob- 
servables. 
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